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ABSTRACT 


The  relaxation  of  a  seeded  plasma  subjected  to  a  constant  electric 
field  is  considered  within  the  formalism  of  the  two-temperature  theory. 
Numerical  solutions  are  obtained  for  streamwise  conductivity  profiles 
in  both  argon  and  nitrogen  plasmas  seeded  with  potassium.  Reaction 
rates  used  here  are  analytical  formulae  based  on  the  Byron  model  for 
collision- induced  recombination.  Extensive  analysis  of  the  computed 
results  and  comparison  with  Zukoski’s  experiments  in  argon  establish 
the  validity  of  the  basic  treatment  and  suggest  the  existence  of  energy 
losses  resulting  from  excitation  of  argon  transitions  originating  at  the 
first  excited  state,  in  addition  to  the  usual  elastic  collision  and  potas¬ 
sium  excitation  energy  losses.  Refined  agreement  between  theory  and 
experiment  requires  consideration  of  corrections  to  the  species- 
averaged  momentum-transfer  cross  section,  after  the  manner  of  Frost; 
however,  present  theories  are  shown  to  overestimate  this  correction 
for  argon.  Details  of  the  Byron  model  are  validated  by  comparison  of 
theoretical  and  experimental  time-resolved  potassium  emission  lines. 
Possible  nonequilibrium  effects  in  nitrogen  are  investigated  with 
various  assumed  molecular  excitation  losses.  Relaxation  times  for 
both  types  of  plasmas  are  correlated  by  use  of  Zukoski's  characteristic 
time  parameter. 
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SECTION  I 
INTRODUCTION 


The  relaxation  of  an  alkali-seeded  plasma  perturbed  by  a  constant 
applied  electric  field  is  of  importance  in  studies  of  nonequilibrium  con¬ 
ducting  plasma  flows  for  two  basic  reasons:  (1)  for  understanding  the 
flow  and  predicting  the  performance  of  magnetohydrodynamic  (MHD) 
energy-conversion  devices  and  (2)  as  a  test  of  the  adequacy  of  the  two- 
temperature  nonequilibrium  plasma  theory  developed  by  various  authors 
and  summarized  in  Ref.  1. 

In  the  first  case,  we  note  that  there  are  conditions  under  which  a 
substantial  portion  of  the  MHD  channel  length  is  required  for  a  per¬ 
turbed  plasma  to  reach  a  steady  state.  These  conditions  not  only  require 
consideration  in  the  design  of  MHD  channels  but  also  influence  possible 
electrothermal  instabilities  (discussed  in  Ref,  1).  Even  when  equilib¬ 
rium  conditions  prevail  in  the  plasma  core,  nonequilibrium  phenomena 
dominate  the  electrode -boundary-layer  region.  In  all  these  MHD  flow 
problems  relaxation  effects  are  of  considerable  interest. 

In  addition  to  its  significance  in  MHD  channel  problems,  the  relaxa¬ 
tion  problem  provides  an  important  test  of  the  two -temperature  theory 
of  nonequilibrium  conduction,  including  the  mechanisms  of  electron- 
energy  loss.  Although  nonequilibrium  electron-fluid  descriptions  abound 
in  the  literature,  Kerrebrock's  (Ref.  1)  assumption  of  a  Maxwell- 
Boltzmann  distribution  of  both  free  and  valence  electrons,  maintained 
at  a  higher  temperature  than  the  bulk  gas  temperature,  reduces  the 
description  of  nonequilibrium  conduction  to  a  tractable  problem  and  pro¬ 
vides  a  means  of  comparing  the  effects  of  various  electron- collision 
energy-loss  mechanisms  with  experimental  conduction  data.  These 
assumptions  are,  generally,  quite  good  for  electron-collision-dominated 
plasmas.  The  experiments  of  Refs.  2  and  3  show  that  satisfactory  agree¬ 
ment  between  experiments  and  theory  can  be  obtained  at  steady-state 
conditions,  provided  proper  account  is  taken  of  the  various  electron- 
collision  effects.  In  the  present  investigation,  detailed  comparison  of 
theoretical  predictions  with  the  relaxation  experiments  reported  in 
Ref.  3  indicates  that  the  two-temperature  theory,  modified  to  include 
ionization  and  recombination,  is  adequate  for  the  full  description  of 
collision- dominated,  relaxing,  nonequilibrium  plasmas. 
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SECTION  II 
BASIC  MODEL 


2.1  ASSUMPTIONS  AND  PROPERTIES  OF  THE  ELECTRON-FLUID  THEORY 


The  electron  fluid  can  be  adequately  described  for  the  present  study 
by  using  conservation  equations  obtained  as  moments  of  Boltzmann’s 
equation,  for  the  electrons  and  their  energy.  For  a  plasma  flowing  at 
constant  velocity  with  applied  electric  field  perpendicular  to  the  stream, 
neglecting  anisotropic  stresses  and  heat  conduction,  the  streamwise 
electron-energy  transport  equation  is 


U 


+  3/2  p 


•) 


je  E  +  I»e 


(1) 


Within  the  restrictions  stated  above  this  equation  is  general,  and 
assumptions  concerning  isotropic  electron-velocity  distribution  affect 
only  the  collision- energy-loss  term  and  the  definition  of  electron 
temperature. 


By  use  of  the  techniques  of  Ref.  4,  one  can  show  that  the  isotropic 
electron-velocity  distribution  is  very  nearly  Maxwellian  at  a  tempera¬ 
ture  higher  than  the  gas  temperature,  provided  electron- electron  col¬ 
lisions  outnumber  collisions  tending  to  distort  the  Maxwellian  speed 
distribution  and  provided  electron  drift  is  much  slower  than  electron- 
thermal  motion.  These  conditions  are  assumed  in  the  present  study. 


The  collision  term  can  now  be  specialized  to  Kerrebrock's  form 
by  use  of  the  expansion  of  Ref.  4. 


-  ~7~  =  3/2  k  (Te  -  T)  ve  St  i»t  Qio  Si 

dt  y  2  yColl  M 


<2) 


dne 

dt 


+  X 


ex  fex  •'ex 


where 


and 


oc 

Si  =  (4/3)  /  e"£  f*  if 


ve 


In  this  equation,  the  first  term  represents  elastic  collision-energy 
loss,  the  second  term  is  the  electron-power  loss  associated  with  ioniza¬ 
tion,  and  the  last  term  corresponds  to  inelastic  collisions  exciting 
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radiative  transitions  in  an  optically  thin  plasma.  Here  S  is  the  conven¬ 
tional  energy-momentum-loss  parameter  -  2  for  Maxwellian  inter¬ 
actions,  8/3  for  hard  sphere  collisions,  and  3.  952  for  electron-ion 
collisions.  The  reference  cross  sections  Qi0  are  Qa  at  0.  4  ev  for 
argon  andQej  for  electron-ion  interactions. 


The  electron- fluid  description  is  completed  by  the  electron- 
conservation  equation: 


In  this  equation  $  is  the  ionization  enhancement  factor  representing  the 
effect  of  multiple  excitation  collisions  as  opposed  to  "one  shot"  ioniza¬ 
tion  collisions.  It  will  be  discussed  in  the  next  section. 


2.2  ELECTRON-COLLISION-DOMINATED  IONIZATION  AND  RECOMBINATION 

A  detailed  discussion  of  ionization  and  recombination  of  hydrogen 
and  of  alkali  atoms  is  presented  in  Ref.  5.  Analytical  formulae  for  the 
reaction  rate  coefficients  are  obtained  by  using  Byron's  assumptions 
(Refs.  6  and  7)  along  with  an  approximate  quantum  mechanical  descrip¬ 
tion  of  the  atomic  structure.  The  results  for  alkali  atoms  agree  with 
numerical  computations  presented  in  Refs.  7,  8,  and  9.  Figure  1 
illustrates  the  basic  properties  of  the  ionization  model  used  in  Ref.  5 
and  in  the  present  study. 

In  this  model  all  states  located  at  or  above  some  atomic  energy 
level  (referred  to  as  the  critical  level)  are  in  Saha  equilibrium  with 
the  free  electrons;  i.  e. ,  their  population  densities  are  governed  by 
the  free  electrons  through  Saha's  equation.  States  at  or  below  one 
principal  quantum  number  below  the  critical  state  are  in  Boltzmann 
equilibrium  with  the  ground  state.  The  rate  of  ionization  or  recom¬ 
bination  is  equal  to  the  net  rate  of  crossing  the  critical  interval.  At, 
defined  by  the  critical  state  and  the  state  one  principal  quantum  num¬ 
ber  below  it. 

•  The  critical  energy  level  is  located  by  minimizing  the  reaction 
rates  with  respect  to  critical  quantum  number.  This  is  done  numeri¬ 
cally  in  Refs.  7,  8,  and  9  and  is  done  analytically  in  Ref,  5  by  using 
approximate  expressions  for  the  average  degeneracy  of  states  within 
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the  critical  interval,  as  a  function  of  critical  quantum  number.  The 
resulting  analytical  values  of  the  critical  energy  level  are  in  good 
agreement  with  Byron's  numerical  values  and  can  be  fairly  well 
approximated  at  moderate  electron  temperatures  by  e  -  fC  =  5k  Te, 
where  <c  is  the  critical  energy  level.  The  critical  interval  is  related 
to  the  critical  energy  level  by 

t^r-  -  2  0  f  (1  +  3/2  #0  (4) 


where 


and 


6 J 


kTe 


C  = 


i  -  fc 

kT. 


The  ionization  enhancement  factor  for  alkali  atoms  is  derived  in  Ref.  5  in 
terms  of  these  parameters: 


3 


(i-s/2  eg 

u TFT 


f  [i  -  »«(i  ^ 


(5) 


where 


and 


ft  = 

1  +N„ 

u  — 

1  + 

r  \ 

0.8 

y.  _u 

e  I  c 

At 

kTe 


An  analogous  formula  for  hydrogen  is  given  in  Ref. 
of  the  atomic  structure  parameters  in  Eq.  (5)  are: 


5. 


Typical  values 


N0  =  1.  85  and  =  1  for  potassium 

N0  =  1.  93  and  =  2  for  cesium 

As  Byron  has  previously  demonstrated  (Refs.  6  and  7),  the  critical 
energy  level  lowers  with  increasing  electron  temperature.  In  the 
formalism  of  Ref.  5,  the  lowest  possible  critical  state  is  the  state  one 
principal  quantum  number  above  ground  state;  therefore,  at  high  elec¬ 
tron  temperatures  Eq.  (4)  must  be  replaced  by 


At  this  limit  G  is  taken  to  be  either  one  or  the  ratio  of  the  average  degen¬ 
eracy  for  all  substates  in  the  lowest  one -quantum  interval  to  that  of  the 
ground  state.  Two  limiting  curves  are  generated  in  this  fashion. 
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In  Fig.  2,  the  computed  three  body  recombination  coefficients  from 
Ref.  5  are  compared  with  other  theories  and  experiments  of  Refs.  9  and 
10.  Equation  (5)  underestimates  the  approach  to  the  limit  for  potassium 
for  reasons  discussed  in  Ref.  5;  however,  the  discrepancy  is  still  well 
within  experimental  scatter.  Equation  (5)  will  be  used  here  with  values  of 
C  from  Ref.  5  for  all  temperatures  below  those  for  which  Eq.  (6)  is  appli¬ 
cable  . 

In  view  of  the  agreement  between  theory  and  experiment  shown  in 
Fig.  2,  it  is  not  the  intent  of  this  report  to  verify  the  reaction  rates  com¬ 
puted  from  the  Byron  model.  Instead,  some  of  the  assumptions  of  the 
model  will  be  examined  critically  by  comparing  experimental  and  theo¬ 
retical  details  of  the  excited- state  relaxation  in  Section  IV. 


2.3  FINAL  RELAXATION  EQUATIONS 


Although  Eqs.  (1),  (2),  (3),  and  (5)  or  (6)  taken  together  with 
appropriate  boundary  conditions  form  a  determinate  set  of  relaxation 
equations,  some  further  simplifications  can  be  introduced  into  the 
present  problem.  The  transport  terms  on  the  left  side  of  Eq.  (1)  are 
orders  of  magnitude  smaller  than  the  terms  on  the  right  side  of  Eq.  (1) 
except  during  an  initial  transient  stated  in  Ref.  3  to  be  about  3  x  10“ 6 
see;  these  terms  will  be  neglected  here.  Physically,  this  means  that 
all  energy  put  into  the  plasma  by  the  electric  field  is  dissipated  by  col¬ 
lisions.  The  electron  temperature  is  determined  by  balancing  energy 
input  against  dissipation.  Additionally,  inelastic  collisions  between 
electrons  and  seed  atoms  will  be  neglected  because  the  findings  of 
Ref.  3  indicate  these  energy  losses  to  be  unimportant  above  current 

densities  of  about  4  am%.  Inelastic  collisions  between  electrons  and 

cm* 

argon  atoms  will  be  absorbed  into  an  effective  energy-loss  collision 
cross  section  in  Eq.  (2).  Equations  (1)  and  (2)  are  now  combined  and 
written  as 


2 

p 


dne 

dx 


+ 


T  “if  k  (Te  '  T)  ni  Si  Qio  ve] 


(7) 


where 


tljkT 

p 


Qi 


The  average  momentum-transfer  cross  section  for  the  "ith" 
is 


Qi(£) 


species 


(8) 


Frost  (Ref.  11)  has  shown  that  errors  are  introduced  by  this  aver¬ 
aging  procedure.  The  effect  of  these  errors  will  be  considered  in 
Section  IV. 
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The  relaxation  is  now  described  by  simultaneous  solution  of 
Eqs.  (7),  (3),  and  (5)  or  {6),  subject  to  the  boundary  condition  that 
the  initial  electron  density  is  some  given  value.  In  practice,  the 
solution  is  obtained  by  finite  difference  integration  of  Eq.  (3)  with 
the  electron  temperature  at  each  value  of  x  being  obtained  by  itera¬ 
tion  of  Eq.  (7)  until  the  specified  electric  field  is  obtained. 


SECTION  III 
COLLISION  INTEGRALS 


The  collision  integrals  appearing  in  Eqs.  (2)  and  (8)  are  of  the 
form  of  Laplace  transforms  of  the  collision  cross  sections.  Unfor¬ 
tunately,  they  can  be  solved  explicitly  only  for  Maxwellian  inter¬ 
actions,  Coulomb  interactions,  and  hard  sphere  collisions.  In  Ref.  3 
these  integrals  were  obtained  numerically  for  each  separate  electron 
temperature.  In  this  report  the  argon  cross  section  integrals  will  be 
computed  by  analytic  integration  of  curve  fits  to  the  cross  sections 
shown  as  functions  of  electron  energy  in  Ref.  3.  Systematic  variation 
of  the  curve  fit  parameters  away  from  the  correct  elastic  collision 
values  permits  simulation  of  inelastic  collisions  exciting  argon  radia¬ 
tive  transitions.  Potassium  and  nitrogen  will  both  be  treated  as  hard 
sphere  gases. 


The  momentum-transfer  integral  can  be  calculated  to  sufficient 
accuracy  by  a  curve  fit  for  the  argon  cross  section  of  the  form 


where 


(9) 


and 

QAo  -  3.82  x  10"‘7  cm1 

The  analogous  cross  section  for  electron-ion  interactions  is  conventional 

Q  =  — £ V  — 

e*  0.582  y  2rrt0  me 

where  A.  is  the  ratio  of  Debye  length  to  the  90 -deg  deflection  impact 
parameter. 
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Provided  bf  0  >  1/2,  the  integrated  argon  cross  section  becomes 

3g  Qaq  (10) 


Qa  = 


8  Ir 


where 


lQ  =  I 


QAo  6 


»-k^0 


Q(f> 


"I’  _(y7  [l  +  erf(j8)] 

vX^+/3)U 

[3/4  +  /3*(3  +  /3s)]  +  /S  (11/2  +  15j8s)| 


and 


0  = 


fafo  -  1/2 


V  f o  n  +  ) 

The  integrated  electron-ion  cross  section  is 


as  used  in  Ref.  3. 


The  energy -loss  integral  in  Eq.  (2)  cannot  be  calculated  from 
Eq.  (9)  because  of  convergence  difficulty.  Instead,  Eq.  (9)  is  used 
for  f  <  f  o»  and  a  parabolic  curve  fit  is  used  for  f  >  f  0-  . 


QA(f) 

Qao 


1  +  a  (f  -  f0)2 


(ID 


The  energy- momentum  parameter,  for  argon  now  becomes 


SA  =  4/3 


-  4 


U  -  b  )’  ” 
1/2 -bfo 


L/2  —  b  f 


1  -  b 


(5/2  +  Z)  +  e 


w 

1  -  b 

(3  +  Zo  +  6Z ) 


1  -  b 


(1  *  Z 


VX  W4  +  +  Zs) 


JtZt  (erf  wz } 

-  erf  (VZ7))  [■+  4/3  [l  +  af0s  T  {3>  -  2af0  r{4  ^  >  +  ar<5i  ]  (12) 


where 

,  (i/2-bfo)J  „  <f0-i/2>3 

*  ( 1  - b ) e0  ’  0  "  d-b)fo 
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and 


r<n-  L ) 


J  e'£  e"_l  A€ 

c 


Matching  the  argon  Ramsauer-Townsend  cross  sections  used  in 
Ref.  3  yields  b£o  =  5.371  and  a£ Q2  =  0.  8752.  The  collision  integrals 
obtained  from  these  parameters  are  shown  in  Figs.  3  and  4  for  com¬ 
parison  with  the  results  of  numerical  integration.  The  analytical 
approximations  are,  obviously,  adequate  for  use  in  the  present  study. 
Also  shown  in  Fig.  4  is  a  curve  obtained  from  b£0  =  5.  371  and 
a£02  =  4.  5.  It  represents  additional  energy  loss  resulting  from  in¬ 
elastic  collisions  and  will  be  discussed  in  Section  IV. 


SECTION  IV 

COMPARISON  OF  THEORETICAL  RESULTS  WITH  ZUKOSKI’S  EXPERIMENTS 
IN  POTASSIUM-SEEDED  ARGON 


The  conductivity  and  current  density  can  be  calculated  from  the 
primary  relaxation  solutions  by  use  of  the  equations  below. 


a 


n  q  kT 

ft  _ 

«n«  ve  Qp 


(13) 


and 


je  m 


(14) 


These  quantities  will  be  compared  with  the  data  reported  in  Ref.  3  in 
the  following  fashion:  First,  the  steady- state  conductivity  will  be 
shown  as  a  function  of  current  density  with  only  elastic  collisions 
accounted  for.  Second,  the  energy-momentum-loss  parameter  will  be 
adjusted  until  theory  and  experiments  agree  well  in  the  steady  state, 
and  the  energy-loss  curve  generated  by  this  procedure  will  be  inter¬ 
preted  quantitatively  as  an  excitation  mechanism.  Third,  the  stream  - 
wise  variation  of  conductivity  with  time  computed  using  this  energy 
loss  will  be  compared  with  the  experimental  relaxation  profiles. 

Fourth,  the  excited- state  behavior  of  potassium  will  be  computed  from 
the  properties  of  the  Byron  model  for  ionization  and  recombination  and 
correlated  with  time -resolved  spectra  reported  in  Ref.  3.  After  the 
adequacy  of  the  basic  treatment  has  been  established,  a  summary  of  the 
relaxation  process  will  be  presented  from  analysis  of  the  relaxation 
solutions. 


As  these  calculations  assume  somewhat  idealized  conditions,  it  is 
well  to  review  the  theoretical  restrictions  and  the  implications  of  effects 
not  taken  into  account  before  comparing  theory  and  experiments. 
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1.  The  theory  assumes  the  electrons  instantaneously  achieve 
equilibrium  with  the  electric  field,  yielding  an  initial  electron  tem¬ 
perature  consistent  with  Eq.  (7)  at  a  given  initial  electron  density. 
Experimentally,  however,  a  finite  time  is  required  for  the  electrons 
to  reach  this  condition. 


2.  Excitation  of  the  seed  has  been  neglected  because  it  was 
shown  in  Ref.  3  to  be  unimportant  at  current  densities  above 

4  g.  This  neglect  causes  some  over  estimation  of  quantities  in 
the  early  relaxation  region. 


3.  Anisotropic  electron -velocity  effects  are  assumed  to  be 
described  adequately  by  the  integral  in  Eq.  (8).  The  rather  drastic 
dependence  of  the  electron- argon  collision  cross  section  on  electron 
energy  should  narrow  the  anisotropic  velocity  distribution  with  respect 
to  that  used  in  Eq.  (8). 

4.  A  related  consideration  is  the  use  of  species -averaged  col¬ 
lision  frequencies  within  the  integral  used  in  Eq.  (8)  instead  of  aver¬ 
aging  over  species  after  performing  the  velocity  integration.  A 
rigorous  consideration  of  these  last  two  effects  is  beyond  the  scope 
of  this  report;  however,  it  will  be  shown  later  in  this  section  that 
present  methods  of  dealing  with  these  two  effects  are  completely 
inadequate  in  the  case  of  seeded  argon.  More  work  of  both  experi¬ 
mental  and  theoretical  character  is  clearly  needed. 


4.1  STEADY-STATE  RESULTS 


Figure  5  shows  that  steady-state  conductivities  computed  with 
only  elastic  collision-energy  losses  (corresponding  to  the  curve  marked 
a£02  =  0-  8752  in  Fig.  4)  agree  only  fairly  well  with  quoted  experimental 
results.  The  present  calculations  for  elastic  collisions  are  exhibited 

/  n  k  \ 

at  only  one  seed  rate  —  =  0.  0023  ]  because  at  all  seed  rates  the  results 

\nA  / 

are  virtually  identical  with  similar  ones  taken  from  Ref.  3.  Better 
agreement  between  theory  and  experiments  is  obtained  by  use  of  the 
energy-loss  curve  (in  Fig.  4)  labeled  af  Q2  =  4.  5.  This  curve  was  gen¬ 
erated  by  forcing  the  computed  conductivity  to  agree  with  experimental 

/nK  \ 

data  at  one  seed  rate  ( —  =  0.  0047  .  Figure  5  shows  that  the  resulting 

\nA  I 

theoretical  conductivities  agree  well  with  experimental  results  at  all 
seed  rates.  A  plausible  mechanism  for  this  additional  energy  loss  will 
now  be  considered  in  detail. 
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Since  the  first  excited  state  of  argon  is  metastable,  its  population 
density  is  higher  than  that  obtained  by  Boltzmann  equilibrium  with  the 
ground  state.  (Its  population  corresponds,  roughly,  to  that  in  the  arc 
column. )  A  likely  electron-energy-loss  mechanism,  then,  is  excita¬ 
tion  of  higher  argon  states  by  collisional  transitions  originating  at  the 
metastable  level.  These  higher  states  should  then  decay  by  allowed 
radiative  transitions.  The  presence  of  metastable  argon  atoms  also 
implies  an  additional  ionization  mechanism  (associated  with  collisional 
de-excitation),  hence  a  mechanism  for  electron-energy  gain.  Including 
both  these  effects,  the  inelastic  collision  energy -momentum-loss 
parameter  for  argon  with  metastable  states  follows. 

(15) 


6a.  ,  .  -  2/3 

■’*  III  C  1  B  StlC 


“e"A7e  dt 


where  « ex  is  the  excitation  potential  for  the  transition  under  considera¬ 
tion,  «  a*  is  the  energy  level  of  the  metastable  state  from  which  the 
transition  originates,  and  n^*  is  the  population  density  of  that  state. 


The  first  term  in  Eq.  (15)  represents  energy  loss  resulting  from 
excitation  of  the  transitions  mentioned.  The  second  term  represents 
energy  put  into  the  electron  gas  by  destruction  of  metastable  states 
resulting  from  argon-potassium  ionization  collisions.  This  term  is 
negligible  except  in  the  early  stages  of  the  relaxation  and  is  omitted 
from  the  relaxation  solutions. 


The  metastable  population  can  be  estimated  by  equating  the  power 
density  supplied  by  the  metastables  to  that  required  to  produce  the 
extrathermal  initial  ionization  observed  in  the  experiments  of  Ref.  3. 
The  initial  conductivities  obtained  in  the  relaxation  experiments  of 

Ref.  3  are  about  0.  5  m  — ,  contrasted  with  gas  thermal  conditions  an 

cm 

order  of  magnitude  lower.  The  electron  temperature  required  to 
produce  such  conduction  is  related  to  the  conductivity  by 


V  8  %/me  kTe  Op 


(16) 


This  corresponds  to  a  required  power  density  obtained  from  .the  steady- 
state  portion  of  Eq.  (7), 


W  =  .  -  -Q  -  ■  pkT  S 

>/ 2trme  kT*  M 


(17) 


The  power  density  supplied  by  the  metastables  is 


VI  =  «* 


(18) 
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where  is  the  cross  section  for  potassium  atom-argon  meta- 

stable  atom  collisions.  It  is  approximately  — Qk;.  The  use  of  these 

equations  yields,  for  a  =  0.48  m  °  and  =  0.  0047,  -A  =  3. 3  x  10"5 

cm  n^  n^ 

corresponding  to  an  arc  column  temperature  of  about  10,  000°K. 


The  only  transition  capable  of  reproducing  the  approximate  shape 
of  the  empirical  inelastic  collision -energy-loss  curve  shown  in  Fig.  6 
(corresponding  to  the  difference  of  the  two  analytical  curves  in  Fig.  4) 
is  IS5  -  2Pg  (Paschen  notation).  Its  excitation  potential  is  1.  623  ev 
(Ref.  12).  When  this  excitation  potential  and  the  calculated  metastable 
population  density  are  substituted  into  the  first  term  of  Eq.  (15),  the 
theoretical  excitation  energy -loss  curve  shown  in  Fig.  6  is  generated. 
In  view  of  the  approximate  nature  of  the  calculation,  the  agreement  of 
these  curves  is  striking.  Future  experimentation  should  certainly  seek 
spectroscopic  evidence  for  these  excitations. 


4.2  CONDUCTIVITY  RELAXATION 

Typical  conductivity -relaxation  profiles,  including  the  effect  of  the 
argon  excitation  energy  losses,  are  shown  in  Fig.  7  in  comparison  with 
experimental  profiles  from  Ref.  3.  The  theoretical  values  have  been 
set  equal  to  the  experimental  values  at  time,  t  =  0,  to  account  for  the 
extrathermal  ionization  present  initially  in  the  experiments.  Also 
shown  are  experimental  and  theoretical  relaxation  times  corresponding 

Q  *”  O'  1 

to  -rrr. — , - „ .  ,  =  1 .  The  agreement  is  satisfactory. 

CTfinal  -  "initial  e 

The  effect  on  the  relaxation  profiles  of  neglecting  excitation  of 
potassium  resonance  radiation  can  now  be  assessed.  In  Ref.  5  it  is 
shown  that  for  cesium  the  ratio  of  radiative  to  collisional  transitions 

1.  24  x  1()14  ,, 

across  the  lowest  lying  critical  interval  is  — - cm-'3.  Since 

ne 

similar  conditions  should  occur  with  potassium,  the  excitation  of 
potassium  resonance  lines  should  be  unimportant  for  electron  densities 
much  greater  than  1.  3  x  1014  cm"3.  The  steady-state  current  density 

corresponding  to  this  electron  density  is  about  2  Zukoski  quoted 

a  current  density  of  about  4  corresponding  to  electron  densities 

of  about  2  x  1014  cm“3.  This  experimental  finding  lends  validity  to 
the  use  of  times  corresponding  to  electron  densities  of  1.  3  x  1014  cm"3 
as  the  criterion  for  separating  regions  of  collisional  and  radiative 
dominance,  as  shown  in  Fig.  7.  Except  for  very  low  electric  fields,  the 
relaxation  profiles  are,  obviously,  collision  dominated,  and  the  present 
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treatment  is  adequate  through  most  of  the  relaxation  region.  The 
observed  discrepancies  between  theoretical  and  experimental  relaxa¬ 
tion  times  can  be  entirely  attributed  to  neglect  of  potassium  reso¬ 
nance  radiation. 


4.3  DETAILS  OF  THE  RELAXATION  MECHANISMS 


Having  established  the  validity  of  the  basic  treatment,  we  can 
now  use  these  solutions  to  shed  light  on  the  excited-state  behavior  of 
the  potassium  seed  and  other  details  of  the  relaxation  process.  Com¬ 
parison  of  theoretical  predictions  with  experimental  excited-state 
populations  obtained  from  time-resolved  spectral  measurements 
reported  in  Ref.  3  provides  a  sensitive  test  of  the  assumptions  asso¬ 
ciated  with  the  Byron  model  for  ionization  and  recombination. 


All  states  located  at  or  above  the  critical  state  are  described  by 

neV  e_lrF7  (19) 


Figure  8  compares  the  spectral  line  intensities  expected  from  states 
described  by  Eq.  (19)  with  experimental  values  from  Ref.  3. 


Looking  first  at  the  resonance  radiation  (1.  61-ev  experimental 
curve),  we  note  that  the  line  intensity  rises  sharply,  reaches  a  plateau, 
and  then  rises  more  slowly  to  its  steady-state  value.  In  contrast,  the 
other  experimental  line  intensity  curve  rises  smoothly  to  its  steady- 
state  value.  There  is  fair  agreement  between  predictions  of  Eq.  (19)  and 
the  observed  radiation  for  the  lower  curve,  indicating  that  the  critical 
energy  level  is  below  3.4  ev. 

The  origin  of  the  plateau  in  the  resonance  radiation  corresponds 
closely  to  the  attainment  of  the  electron  density  (1.  3  x  10 14  cm-3)  at 
which  collisional  and  radiative  transitions  are  equal  across  the  lowest 
one-quantum  interval  above  ground  state.  The  occurrence  of  this 
plateau  and  the  line  intensity  behavior  after  the  plateau  suggests 
strongly  that  the  first  excited  state  of  potassium  (located  at  1,61  ev) 
lies  within  the  critical  interval. 


These  considerations  restrict  the  critical  energy  level  at 
Te  =  3100°K  to  the  range  from  2.  61  to  3.06  ev.  Reference  5  predicts 
«c  =  2.  67  ev;  the  nearest  state  satisfying  the  assumption  of  a  one- 
quantum  critical  interval  is  5S,  lying  at  2.61  ev.  References  5  and  8 
for  cesium  imply  that  the  critical  state  in  potassium  should  be  5P, 
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corresponding  to  * c  =  3.06  ev.  The  Byron  model  is  fairly  well  sub¬ 
stantiated  in  spite  of  the  dubious  nature  of  some  of  the  statistical 
assumptions  at  high  electron  temperature. 


The  dip  in  the  electron  temperature  shown  in  Fig.  9  can  be 
explained  on  the  basis  of  competition  between  ionization  and  other  col¬ 
lision  phenomena.  Let  the  ratio  of  power  loss  by  ionization  to  that 
resulting  from  all  other  energy-loss  collisions  be 


S 


£ 


3/2  -^7e  k(T 

M 


e 


T>S|»,Qi„Si 


Equation  (7)  can  be  rewritten  as 


L  2  _me.  T  e  / T  „  A  /  p  (A 
n  M  T  \  T  "/  \  <!e  / 


8  (1  +  S) 


(20) 


where 


8  = 


Sini  QisSi 

IFF  Q 


The  most  rapidly  varying  quantities  in  Eq.  (20)  are  S_and  Q;  S  decreases 
with  time  as  recombination  becomes  important,  and  Q  increases  with 
time  as  Coulomb  scattering  occurs.  For  an  ionization-dominated  plasma, 
Eq.  (20)  predicts  that^the  electron  temperature  should  rise  monotonically 
until  the  increase  of  Q  causes  it  to  level  off.  In  contrast,  a  plasma  whose 
ionization  losses  are  much  smaller  than  other  losses  (such  as  large  ex¬ 
citation  losses)  should  exhibit  an  electron-temperature  decline  to  steady- 
state  conditions. 


The  relaxation  of  seeded  argon  exhibits  both  these  types  of  behavior. 
Figure  9  shows  that  the  relative  electron-ion  collision  frequency  in¬ 
creases  faster  than  S  decreases,  causing  the  electron  temperature  to 
decrease  until  the  two  terms  balance.  Beyond  this  time,  recombination 
proceeds  faster  than  the  increase  of  electron-ion  scattering,  causing  the 
electron  temperature  to  increase  to  its  steady-state  value.  This  behavior 
is  typical  of  plasmas  in  which  ionization  energy-loss  rates  are  comparable 
to  all  other  energy-loss  rates. 


4.4  UNCERTAINTIES  ASSOCIATED  WITH  COLLISION  CROSS  SECTION  AVERAGING  PROCEDURE 


The  preceding  discussion  has  shown  that  moderate  density,  alkali- 
seeded  plasmas  are  well  described  by  the  two-temperature  theory.  In 
particular,  the  assumption  of  a  Maxwell- Boltzmann  isotropic  distribution 
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of  free  electrons  and  bound-state  electrons  (above  the  critical  state 
during  relaxation)  is  adequate  for  the  description  of  most  collision- 
dependent  phenomena.  Unfortunately,  effects  dependent  on  the 
anisotropic  velocity  distribution  can  not,  at  present,  be  predicted 
with  such  certainty. 


The  averaged  collision  cross  section  in  Eq.  (7)  should  really  be 
obtained  from 


~  =  4/3  V" 
PQ 


€  fe*  (£)  i€ 


(21) 


When  the  approximate  solution  of  Boltzmann1  s  equation  given  in  Ref.  4 
is  good,  the  ratio  of  anisotropic  to  isotropic  distributions  is 


-qe  E  dii.fo°(£)  l 

=  kT.  d£ 


(22) 


and  Eq.  (21)  can  be  rewritten  as 


ni  Qi(f) 


where 


H(£> 


dfe°  H)  /  d(e°  (£) 


Maxwellian 


(23) 


When  scattering  by  any  one  species  predominates,  Eq.  (23)  reduces  to 
the  cross  section  used  in  Eq.  (7),  the  form  first  assumed  in  Ref.  13  to 
correlate  shock  tube  conductivity  studies  in  unseeded  argon. 


For  the  case  of  H  =  1,  Frost  (Ref.  11)  has  shown  that  when  scatter¬ 
ing  by  ions  is  comparable  to  neutral  scattering,  the  conductivity  is  over¬ 
estimated  by  use  of  the  cross  sections  in  Eq.  (8)  and  must  be  ratioed 
down  by  use  of  a  correction  factor.  Figure  10  shows  the  conductivity- 
correction  factor  (ratio  of  cross  section  given  by  Eqs.  (8)  and  (23)) 
obtained  from  Eq.  (23)  (with  H  =  1)  for  a  plasma  composed  of  hard 
spheres  and  ions.  It  is  slightly  lower  than  Frost's  because  Frost  used 
a  different  ion  cross  section  (corresponding  to  a  different  form  of  the 
electron -electron  scattering  term);  however,  the  curves  are  similar. 
Also  shown  are  results  for  nitrogen  obtained  by  using  curve  fits  of 
experimental  cross  sections  in  an  analytic  calculation. 


When  Eq.  (23)  was  evaluated  numerically  for  partially  ionized, 
seeded  argon,  the  corrections  were  found  to  be  more  severe  than  those 
mentioned  above  -  so  severe  that  the  conductivities  shown  in  Fig.  5 
were  underestimated  by  almost  a  factor  of  2.  As  this  result  is  incon¬ 
sistent  with  the  good  correlation  of  experimental  data  reported  in 
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Refs.  2,  3,  and  13,  we  are  led  to  conclude  that  Eq.  {23)  is  inadequate 
in  argon.  Two  possible  reasons  for  this  inadequacy  will  now  be 
explored. 

First,  consider  the  effect  of  electron-argon  collisions  on  the 
electron-distribution  function:  The  drastic  dependence  of  the  argon 
collision  cross  section  on  electron  energy  {resulting  from  the 
Ramsauer-Townsend  effect)  narrows  the  anisotropic  velocity  dis¬ 
tribution.  If  the  limit  form  H(£)  e“£  — *-6(f  -  £  o)  used  as  the 
distribution  function  in  both  Eqs.  (23)  and  (8),  the  two  equations  pre¬ 
dict  the  same  cross  section,  suggesting  that  proper  results  of  using 
Eq.  (23),  if  the  velocity  distribution  were  known  correctly,  should 
be  closer  to  those  of  Eq.  (8)  than  of  Eq.  (23)  with  the  usual  assump¬ 
tion  that  H  =  1. 

An  alternate,  but  related,  possibility  is  the  following*:  Occasional 
electron-argon  collisions  in  a  sequence  of  deflections  resulting  from 
electron-ion  interaction  may  drastically  change  the  scattered  velocity 
distribution,  affecting  subsequent  electron- ion  collisions .  Scattering 
under  such  conditions  is  similar  to  scattering  resulting  from  a  com¬ 
bination  of  Coulomb  interactions  and  inelastic  collisions.  In  such 
circumstances  it  is  probably  appropriate  to  perform  velocity- 
distribution  averages  before  averaging  over  species,  leading  to  the 
use  of  Eq.  (8)  instead  of  Eq.  (23).  Rigorous  determination  of  the 
factor  H  likely  will  require  a  numerical  solution  of  the  Boltzmann- 
Fokker- Planck  equation,  similar  to  that  mentioned  in  Ref.  11  (for  a 
gas  composed  of  hard  spheres,  electrons,  and  ions).  It  is  questionable 
whether  such  sophistication  at  this  stage  is  worthwhile. 


SECTION  V 

THEORETICAL  RESULTS  FOR  SEEDED  NITROGEN 


Because  of  its  importance  in  practical  MHD  devices,  relaxing  con¬ 
duction  in  seeded  nitrogen  will  now  be  studied.  Lack  of  reliable  experi¬ 
mental  conductivity  data  and  lack  of  understanding  of  energy-loss 
mechanisms  will  necessarily  restrict  the  scope  of  this  section.  It  is 
intended,  however,  that  this  study  should  guide  future  experimentation 
in  seeded  diatomic  plasmas.  With  this  purpose  in  mind  we  shall  study 
relaxation  in  nitrogen  at  only  one  seed  mole  fraction  but  at  various  gas 
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temperatures.  Calculations  will  be  presented  using  two  forms  of 
electron- energy  loss  associated  with  molecular -state  excitation. 

Electronic  collisionai  excitation  of  diatomic  molecules  is  not 
well  understood.  Lacking  reliable  excitation  cross  sections,  we 
shall  lump  electron-neutral  energy-loss  collisions  and  arbitrarily 
use  values  of  25  and  100  for  the  resulting  parameter 

8/3  nR  Qk  +  Sy  nN  Qn 

°o  -  - Z - 

p  Q 
ItT 

The  first  value  corresponds  to  experimental  nitrogen  excitation 
losses  at  room  temperature  (Ref.  14).  The  second  value  represents 
a  crude  attempt  to  account  for  additional  excitation  at  the  high  gas 
temperature  characteristic  of  current  MHD  devices.  Experimentally, 
50  =  100  correlates  the  data  of  Ref.  15;  however,  uncertainties  in  the 
gas  temperature  measurement  and  electrode  sheath  and  boundary - 
layer  voltage  drops  obscure  any  conclusive  analysis  of  the  results  of 
Ref.  15.  It  is  probably  best  to  view  60  =  100  as  an  estimated  upper 
limit  for  the  electron-energy  loss  in  seeded  nitrogen. 

Some  theoretical  justification  for  this  value  can  be  obtained  by 
treating  electron- collisionai  excitation  of  an  anharmonic  oscillator. 
Because  upper  vibrational  energy  levels  are  more  closely  spaced 
than  lower  ones,  more  one-quantum  excitations  are  possible  at  high 
gas  temperature  than  was  true  in  the  room  temperature  experiments 
reported  in  Ref.  15.  Assuming  classical  collisionai -excitation 
probability  and  using  energy  levels  for  the  anharmonic  oscillator 
(Ref.  16),  one  can  show  that  the  inelastic  portion  of  is  approxi¬ 
mately  proportional  to  -\J t7  This  order  of  magnitude  estimate  pre¬ 
dicts  60  =•  75. 

Figure  11  shows  both  steady -state  and  relaxing  conductivity  cal¬ 
culations.  It  is  immediately  apparent  that  experiments  in  an  unrelaxed 
flow  can  lead  to  erroneously  high  apparent  values  of  60,  if  the  relaxa¬ 
tion  is  assumed  complete.  It  should  also  be  noted  that  experimental 
determination  of  60  from  the  curvature  of  the  conductivity-current 
density  curves  will  require  very  precise  conductivity  measurement  - 
likely  10  percent  or  better.  Probably,  spectroscopic  investigations 
will  be  required  to  determine  50  accurately. 

Although  they  are  not  shown  here,  the  relaxation  profiles  in  seeded 
nitrogen  are  very  similar  to  those  in  seeded  argon.  In  nitrogen,  how¬ 
ever,  the  electron  temperature  always  stays  constant  or  monotonically 
decreases  with  time  because  the  large  molecular  excitation  losses 
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dominate  the  relaxation  process,  as  indicated  in  Section  4.  3.  This 
relaxation  behavior  is  in  marked  contrast  to  that  noted  in  argon  because 
of  the  difference  in  collisional- energy -loss  magnitudes. 

The  relaxation  times  in  both  nitrogen  and  argon  are  shown  in 
Fig.  12  as  functions  of  Zukoski's  characteristic  time  parameter, 

.  ^‘initial  f 

*c  =  T5 - 

^initial  E 

The  nonlinear  curves  for  nitrogen  indicate  that  relaxation  in  nitrogen 
is  dominated  by  the  large  energy  losses  in  contrast  to  the  nearly  linear 
curves  for  argon,  where  the  relaxation  is  power  limited.  Since  flow 
residence  times  in  many  MHD  channels  are  of  the  order  of  200  micro¬ 
seconds,  the  relaxation  times  in  nitrogen  at  temperatures  below  2500<>K 
seem  prohibitively  large. 


SECTION  VI 
CONCLUDING  REMARKS 


This  report  has  considered  the  relaxation  of  a  moderate  density, 
alkali-seeded  plasma  subjected  to  a  constant  electric  field.  By  corre¬ 
lating  numerical  solutions  with  experimental  data  in  argon  reported  in 
Ref.  3,  the  adequacy  of  Kerrebrock's  two -temperature  electron-fluid 
theory  in  combination  with  Byron's  ionization  and  recombination  model 
has  been  established.  In  particular,  fair  agreement  between  theoretical 
and  experimental  conductivity  is  obtained  when  only  elastic  collision 
energy  loss  is  considered.  Allowance  for  inelastic  collisions  improves 
the  agreement.  Detailed  analysis  of  time -resolved  spectra  from 
potassium -excited  states  indicates  that  the  excited-state  behavior 
assumed  in  the  Byron  model  is  approximately  correct  and  that  the 
critical  state  is  correctly  located. 

Correlation  of  experimental  initial  conditions  of  Ref.  3  (extra- 
thermal  ionization)  with  the  presence  of  argon  metastable  atoms  indi¬ 
cates  that  a  likely  inelastic  collision  electron- energy-loss  mechanism 
is  the  excitation  of  transitions  originating  at  the  metastable  first 
excited  state.  Future  experimentation  in  argon  should  include  a  spec¬ 
troscopic  verification  of  this  mechanism. 

The  procedure  for  obtaining  averaged  momentum  transfer  collision 
cross  sections  is  questioned.  It  is  shown  that  better  results  are  obtained 
in  argon  by  use  of  the  simple  collision-frequency  addition  rule  of  Ref.  13 
than  by  more  sophisticated  calculations  similar  to  those  of  Ref.  11.  The 
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latter  likely  fail  because  of  inconsistencies  in  computing  the  aniso¬ 
tropic  velocity  distribution  for  the  case  in  which  an  electron's  energy 
changes  appreciably  during  deflection  through  90  deg. 

Relaxing  and  steady-state  conductivities  in  seeded  nitrogen  have 
been  computed  for  two  values  of  electron- energy  loss  caused  by 
excitation  of  molecular  energy  levels.  It  is  shown  that  very  precise 
experiments  will  be  required  to  establish  conclusively  these  losses 
from  conductivity  measurements  alone.  It  is  suggested  that  spec¬ 
troscopic  studies  be  done  in  this  area,  also.  Comparison  of  relaxa¬ 
tion  times  in  nitrogen  and  argon  shows  argon  to  be  power  limited, 
whereas  the  rate  of  relaxation  in  nitrogen  is  strongly  affected  by  the 
large  molecular  excitation  losses. 
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Fig.  1  Model  of  Alkali  Atomic  Structure  Used  in  Reaction  Rate  Calculations 
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Fig.  4  Energy-Momentum-Loss  Parameters  for  Argon 
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Fig.  6  Inelastic  Collision  EnergyMomentum-Loss  Parameters  for  Argo 
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Fig.  8  Theoretical  and  Experimental  Potassium  Excited-State  Relaxation  Profiles 
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Fig.  10  Corrections  to  Conductivity  According  to  Eq.  (23),  (H  *=  1 ) 
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Fig.  12  Correlation  of  Relaxation  Time  with  Characteristic  Time  Parameter 
for  Power-Limited  Relaxation 
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